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ABSTRACT 

We present near-infrared (NIR) adaptive optics imaging obtained with 
VLT/NACO and optical spectroscopy from the Southern African Large Telescope 
(SALT) of a luminous infrared galaxy (LIRG) IRAS 19115-2124. These data are com- 
bined with archival HST imaging and Spitzer imaging and spectroscopy, allowing us 
to study this disturbed interacting/merging galaxy, dubbed the Bird, in extraordinary 
detail. In particular, the data reveal a triple system where the LIRG phenomenon is 
dominated by the smallest of the components. 

One nucleus is a regular barred spiral with significant rotation, while another is 
highly disturbed with a surface brightness distribution intermediate to that of disk 
and bulge systems, and hints of remaining arm/bar structure. We derive dynamical 
masses in the range 3-7xlO lo M0 for both. The third component appears to be a 
1-2x10 10 Mq mass irregular galaxy. The total system exhibits HII galaxy-like optical 
line ratios and strengths, and no evidence for AGN activity is found from optical or 
mid-infrared data. The star formation rate is estimated to be ~ 190 M^yr -1 . We also 
report a search for supernovae from NIR images separated by 5 months and search 
for super star cluster candidates. We detect outflowing gas from the Bird mostly in 
the range 100-300 km s _1 using Nal D absorption features. Overall, the Bird shows 
kinematic, dynamical, and emission line properties typical for cool ultra luminous IR 
galaxies. However, the interesting features setting it apart for future studies are its 
triple merger nature, and the location of its star formation peak - the strongest star 
formation, as revealed by Spitzer imaging, does not come from the two major A-band 
nuclei, but from the third irregular component. This is in contrast to the conventional 
view that the (U)LIRG phases are powered by infalling gas to the major nuclei of the 
merging spiral galaxies. Aided by simulations, we discuss scenarios where the irregular 
component is on its first high-speed encounter with the more massive components. 

Key words: galaxies:individual(IRAS 19115-2124) - infrared: galaxies - galaxies: 
starburst - galaxies: interactions - galaxies: evolution - galaxies: kinematics and dy- 
namics 
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1 INTRODUCTION 

Over the last three decades, evidence has mounted that lu- 
minous infrared (IR) galaxies (LIRGfl) ma y signpost im- 
portant events in understanding the wider picture of galaxy 
formation and evoluti on. First hinted at by ground based 
studies in the 1970s l|Rieke fc Lowlll972l ). the importance 
of infrared-luminous activity in galaxies was clearly demon- 
strated by surveys with IRAS (frlouck et al.l ll985T l. which 
found large numbers of extragalactic sources that were 
bright at mid/far-IR wavelengths, and varied substantially 
in nature as a function of infrared (8-1000/im) luminosity. 
Those sources with IR luminosities less than about 1O 1O L0 
were mainly either dwarf galaxies or relatively dustless ellip- 
ticals, where the IR emission arose almost entirely from dust 
grains heated by ambient interstellar light. Sources with IR 
luminosities in the range 10 10 L Q < L ir < 10 115 L Q con- 
tained a significant fraction of ellipticals, but with a rising 
number of IR-luminous disk galaxies and galaxy mergers. 
At IR luminosities of > 10 11,5 Lq however, mergers were 
almost ubiquitous, with very few undisturbed disk or el- 
liptical systems. In the bulk of these systems the IR emis- 
sion was thought to arise from some combination of dense, 
compact starbursts, and d eeply buried AGN. Surveys with 
ISO ( e.g. lElbaz et alJliool ) and SCUBA (e.g. lHughes et all 
Il998l : ICoppin et al.ll2006h then showed that LIRGs become 
increasingly common with increasing redshift, going from 
of order a few hundred examples over the whole sky at 
z < 0.1 to several hundred per square degree at z ^ 1. 
These distant LIRGs seem similar, a t least superficially, to 
their low redshift counterparts (e.g. iFarrah et aLll2002alll 
Chapman et al.ll2003l : ISmail et alj|2003l ; ' 



Takataet al. 200 



Valiante et alj 2007 ). Excellent reviews of the pr operties of 



LIRGs can be found inlSanders fc Mirabe and more 

recently in lLonsdale et all l|200rj) . 

That the luminous end of LIRGs were discovered to 
be almost invariably mergers containing obscured starbursts 
and AGN linked them to several important galaxy trans- 
formational processes, but also raised many questions. Nu- 
merical simulations suggested that mergers could s erve to 
trans form disk galaxies into elliptical galaxies (e.g. iBarned 
1990), potentially linking LIRG activity to a key step in 
the formation of the familiar galaxy morphologies seen lo- 
cally, but some recent simulations have found evidence 
that this may not always be the case (e.g. iBournaud et alj 
120051 ; [Robertson et alj|2006l ). The exact number of merger 
progenitors is also controversial, with some studies find- 
ing that LIRGs are mostly mergers between two galaxies 
l|Veilleux et al|[2002l ). and others finding eviden ce for multi- 
ple mergers in a signifi cant number of LIRGs (|Borne et al.l 
l200d ; lFarrah et aj|200l1 ). The starbursts and AGN in LIRGs 
are obviously candidates for building up the stellar masses 
and large central black holes seen in ellipticals, but ex- 
actly how and when starburst and AGN activity are trig- 
gered as a function of merger stage, how the star formation 



is distributed spatially and evolves in time, and how the 
starburst and AGN contributions to the total IR emission 
can be disentangled, are still debated. Two examples in- 
clude the l ongstanding debate ov er the links between LIRGs 
and QSOs (Kawakatu et al. 200^), and the recent suggestion 
that the far-IR emission from LIRGs is determined more 
by whether the starburst or AGN is the dominant feedback 
mechanism, rather than w hich is more intrinsically luminous 
IjChakrabarti et alj 12007). The issue of feedback in LIRGs, 
in particular large-scale outflows, is also important for de- 
termining how LIRGs ma y contribute to the enrichment of 
the in tergalactic medium (|Veilleux et alj [20031 ; iLfpari et alj 
120051) . 

The central position that LIRGs seem to play in galaxy 
evolutionary processes at all redshifts makes their in-depth 
understanding of great interest. When considering studies of 
local LIRGs there are, broadly, two approaches; population 
studies of large samples to assess trends, and concentrated 
studies of single objects using panchromatic datasets. In this 
paper we adopt the latter approach, utilizing a comprehen- 
sive library of imaging and spectroscopic data to study the 
dynamics and power source within a nearby gas-rich inter- 
acting LIRG, to understand how merger dynamics are affect- 
ing, and have affected, its starburst and AGN activity. Our 
dataset comprises imaging from the Advanced Camera for 
Surveys onboard the Hubble Space Telescope (ACS/HST), 
deep optical spectroscopy from the newly commissioned 
Southern African Large Telescope (SALT), adaptive optics 
near-infrared (NIR) imaging from the Very Large Telescope 
( VLT) , and mid-infrared photometry and spectroscopy from 
the instruments onboard Spitzer. Using this dataset we (1) 
construct a detailed picture of the merger dynamics in this 
system, (2) estimate the number of merger progenitors, (3) 
deduce the luminosities and other relevant parameters for 
starburst and AGN activity, (4) search for evidence for out- 
flows, and (5) discuss scenarios of the system's evolution. 
We assume a spatially flat cosmology with Ho = 73 km s~ x 
Mpc" 1 , Q, = 1, and Q m = 0.27. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 The target 

The luminous infrared galaxy IRAS 19115-2124, also known 
as ESO 593-IG 008, lies at a redshift z = 0.04 9. It has an 
IR luminosity of 10 119 L© (ISanders et alj|2003l ). making it 
fall just short of a formal ULIRG classificati on of 10 12 Lp), 
and a molecular hydrogen mass of3xlO lo M l|Mirabel et all 
1990). It is a classical example of a nearby, IR-luminous, gas- 
rich interacting system. The Virgo infall corrected distance 
to IRAS 19115-2124 with our adopted systemic velocity of 
14576 ± 9 km s _1 is 200.0 Mpc, and the angular scale 0.97 
kpc arcsec - . 



1 Based on observations made with ESO Telescopes at the 
Paranal Observatory under programme 073.D-0406A, and with 
the Southern African Large Telescope (SALT). 

2 objects with 8— lOOO^tm luminosities in excess of 10 11 Lq; ob- 
jects with 8-1000^tm luminosities > 10 12 Lq are commonly re- 
ferred to as ultraluminous IR galaxies, ULIRGs 



2.2 NACO near-infrared imaging 

Near-infrared imaging of IRAS 19115-2124 was carried out 
using the NAOS-CONICA (NACO) adaptive optics (AO) 
instrument on the VLT UT4. The Ks-band (the subscript 
will be dropped henceforth for simplicity) data set, using 
the S27 camera giving a field of view (FOV) of 27 arcsec 
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and pixel size of 0.027 arcsec, presented here was taken in 
service mode on April 13, 2004. The visual wavefront sen- 
sor in mode VIS-WFS 2-7 using a V = 12.9 mag natural 
guide star 22 arcsec southwest of the target was used for the 
AO correction. Coherent energies of the reference star in the 
range 48 to 54 per cent were achieved in individual frames, 
while a strehl ratio (SR) ~ 0.25 was measured at the loca- 
tion of the target galaxy in the final image. The FWHM of 
point sources close to the target galaxy in the final image are 
approximately 0.10 arcsec. The NACO if-band observations 
were repeated on September, 10, 2004, to allow searching for 
supernovae (SNe) within the nuclear regions of the galaxy. 
This time coherent energies for the reference star ranged be- 
tween 36 and 51 per cent in the individual frames, yielding 
SR~0.18 at the location of the target galaxy in the final 
image. 

The (auto)jittered frames were median combined to 
form a sky image using IRAF. The sky subtracted images 
were then de-dithered using the centroid of a bright field 
star. The AO correction was sufficiently stable over the 
whole observation that 23 out of the 24 frames obtained 
on April 13 were used to median combine the final image. 
The individual frames had 90 sec exposure times, giving a 
total integration time of 2070 seconds for the final image of 
IRAS 19115-2124 (Fig. [[J. For the second epoch image all 
24 of the 90 sec exposures taken were median combined to 
form the final frame with 2160 sec total integration time. 

We tied the photometry of the final image to the 2MASS 
if-band using large-aperture photometry of the galaxy itself, 
as well as 2MASS stars in the FOV. The resulting zero-point 
is also consistent with the ESO provided photometric stan- 
dard star calibration, and we estimate the absolute Vega- 
based magnitude to be accurate to within ±0.07 mag. The 
astrometric solution was found using GSC2 field stars. Four 
were available in the FOV, and the rms of the fit was 0.09 
arcsec; the absolute GSC2 system has an uncertainty in the 
range 0.2-0.5 arcsec for individual stars. 

2.3 Optical observations 

2.3.1 HST/ACS imaging 

We extracted archival ACS images (PI: Evans) of the Bird 
in F435W and F814W bands, i.e. in the HST B and I- 
bands, of exposure times 1260 and 720 seconds, respectively. 
These data are already calibrated and we used the supplied 
Vega-based zeropoints. We did not perform any additional 
processing of the images, except a flux-conserving geometric 
transformation into the NACO image size and resolution for 
the purpose of photometry in identical apertures. The im- 
ages, as well as a 3-colour image with the if-band included, 
are shown in Fig. [T] 

2.3.2 SALT/RSS spectroscopy 

Spectroscopic observations presented here w ere obtained 
with the Robert Stobie Sp ectrograph (RSS, iBurgh et al.l 
2003; lKobulnickv et 811 2003) during th e commissioning and 
performance verification stage of SALT (jBucklev et al . 2006; 
lO'Donoghue et al.l [2006), a new 10-m class telescope in 
Sutherland, South Africa. 

The RSS observations are described in Table [1] We 



Table 1. Log of VLT NIR imaging and SALT optical spec- 
troscopy of IRAS 19115-2124. 





Date 


Exp. time 


Band 




VLT/NACO 


13/04/2004 


2070 sec 


K s 




VLT/NACO 


10/09/2004 


2160 sec 


Ks 






Date 


Exp. time 


Slit 


PA 


SALT/RSS 


19/09/2006 


1500 sec 


1.5" 


192 deg 


SALT/RSS 


10/10/2006 


900 sec 


1.0" 


108 deg 


SALT/RSS 


22/10/2006 


1200 sec 


1.5" 


90 deg 



placed slits through the galaxy system in three different 
position angles (PA) on three different nights. The slit po- 
sitions are indicated in Fig[T] The same grating (PG1800) 
with a dispersion of 0.41 A pix -1 and spectral resolution of 
2.4 A FWHM (with 1.5 arcsec slit), and spectral range from 
« 5860 A to » 7150 A was used each time. The October 22 
data set is of the best quality, taken in both good seeing and 
transparency conditions. The October 10 observations were 
taken through thin clouds, and have lower signal to noise 
ratio (SNR) than the others, though all the strong emis- 
sion lines are very clearly seen. Spectra of Cu-Ar and Ne 
comparison lamps were taken after the science frames. 

The data were first bias and overscan subtracted, gain 
corrected, trimmed, cross-talk corrected, and mosaiced, us- 
ing a specially adapted salt IRAF package, and cosmic ray 
corrected using MIDAS routines. The IRAF twodspec pack- 
age was used to wavelength calibrate and correct each frame 
for distortion and tilt, and to background subtract the 2D 
spectra. The apall routine was used to extract various ID 
spectra. The derived internal errors for the 2D wavelength 
calibrations are a — 0.05 A throughout the wavelength range 
of detected emission lines, i.e. < 2.5 km s^ 1 , except for 
September 19 where a — 0.09 A. In addition, a check of night 
sky lines thr oughout the wavel ength range, using methods 
described in IZasov et all (|200Ch , confirmed the accuracy of 
the calibration to this level. All line-of-sight velocities pre- 
sented here were then corrected for heliocentric motion. 

Figure [2] shows the ID and 2D spectra of the PA=90 
data, while Fig. [3] shows the 2D spectral region around Hct 
of all three slit positions. 



2.4 Spitzer MIR imaging and spectroscopy 

We obtained archived Spitzer IRAC (3.6, 4.5, 5.8, and 8um) 
and MIPS (24^m) images of IRAS 19115-2124 (PI: Maz- 
zarella). We used the post-basic calibrated data (PBCD) 
products provided by the S14.0.0 and S14.4.0 versions of 
the pipeline for the IRAC and MIPS data, respectively. The 
IRAC PBCD mosaics were rebinned to the pixel scale of 
NACO and were aligned to the same orientation and pixel 
coordinates using centroid positions of 2-4 point sources vis- 
ible in both images. Fig. [3] shows the IRAC and MIPS con- 
tours. 

Spitzer IRS data (PI: Armus) were also extracted from 
the archive. The data consists of SL (5.2 - 14.5 urn wave- 
length range, 3.7 arcsec slit width) and LL-mode (14 - 38 
urn, 10.5 arcsec slit width). Figure [5] shows the spectrum 
with the most important polycyclic aromatic hydrocarbon 
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Figure 1. The NACO image of IRAS 19115-2124 (obtained on April 13, 2004) is at the top, and the three slit positions observed with 
SALT/RSS are indicated, as are the main components of the interacting system, discussed and named in Section 13,1.11 At lower left 
are the HST B and I images, and at lower right the combined BIK 3-colour image. All tick marks are in one arcsec intervals, and the 
brightness scales are logarithmic, except for the 3-colour image. 
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Figure 2. The reduced 2D spectrum around the target at PA=90 covering nearly the whole observed wavelength range (cut slightly 
in the blue). The ID spectrum is extracted from within a wide 10 arcsec aperture. The two gaps in the spectrum are gaps in between 
the 3 CCD chips of the RSS detector, and the high noise area at w 6720 A is an incomplete subtraction of scattered light from SALT 
'autocollimator laser'. 



(PAH) and emission lines indicated. The orientations of the 
SL and LL mode slits are overlaid on the IRAC and MIPS 
contour images, respectively, in Fig. [4] 



The IRS data were reduced using the SMART and 
SPICE software packages. The individual BCD frames for 
each nod position were cleaned of cosmic rays and other ar- 
tifacts using the irsclean software package, and combined 
into a single image using the 'fair coadd' option within 
SMART. Sky subtraction was performed for each module by 
using the off-order image from the same module (so called 
'order-order' sky subtraction). Spectra were then extracted 
and flux-calibrated from these sky-subtracted, combined im- 
ages using the 'optimal extraction' option within SPICE. 
All other parameters, such as extraction apertures (> 14 
arcsec), were left at their default values. The ID spectra for 
each nod position were then combined to give the final spec- 
trum. However, since the IRS slits do not completely cover 
the galaxy system, and because of the varying slit size in 
the short and long wavelength part of the IRS spectrum, 
we multiplied the SL section of the this final spectrum by 
a factor of 1.6 and the LL section by 1.2. These corrections 
were estimated from flux falling outside the slits in the IRAC 
and MIPS maps, and we estimate them to be accurate to 
~ 10-15 per cent. The IRS spectrum agrees with broad band 
photometry after the correction. 



3 ANALYSIS 

3.1 System characteristics from imaging 

3.1.1 General morphology and photometry 

On the basis of the if-band NACO AO-image alone, IRAS 
19115-2124 appears to be a major merger system with at 
least two prominent components, perhaps three, with tidal 
tails in addition (see Fig. [TJ. With its outspread wings, a 
head complete with a beak, a body and a tail, we dub this 
merging galaxy system the Bird. 

We present photometry of the major components in 
the optical and NIR b ands below. For these SExtractor 
( Bcrt in fc Arnoutd 119961 ) and IDL based aperture photom- 
etry was used, and all the magnitudes are Vega based. We 
note that no aperture corrections, or corrections dependent 
on the slightly varying SR ratio over the image are needed 
when the aperture used is greater than 0.5 arcsec in radius, 
or in the case of extended objects. Photometry of smaller 
point-like sources is explained in Sections I3.1.4l and l3.1.3l 

The largest and brightest component in Jf-band is the 
Southern galaxy nucleus, or the Body of the Bird. It is 
severely extinguished and virtually invisible in the ACS B- 
band image. The bright NIR extent of this galaxy is ~ 2 kpc. 
In a 1 kpc radius aperture (1.03 arcsec) the apparent bright- 
ness of the central region is 7^=12.82 mag, corresponding to 
Mk = —23.69 absolute magnitude. This is not much fainter 
than an L* galaxy, M% = -24.1 ICole et all l|200ll ). Mag- 
nitudes in other apertures and all bands are tabulated in 
Table [2] The galaxy is quite disturbed and extends South 
with a more than 6 kpc long Tail of the Bird, which has 
several clear non-pointlike, elongated 0.3-1.0 kpc sized sub- 
concentrations. The brightest two of these have K ~ 17.0 
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Figure 4. Spitzer IRAC 3.6, 4.5, 5.8, and 8.0 um, and MIPS 24 \xm contours are overlaid on the NACO _fC-band data. Contour levels 
are logarithmic. The peak flux from Spitzer bands shifts from the Body towards the Head NIR component as the wavelength increases. 
The dotted lines on IRAC panels denote the width and orientation of the slit used in IRS LS spectrum (5.2-14.5 um), and the MIPS 
panel is overplotted with the LL slit (14-38 um). 



(B = 20.9, / = 20.0 mag), i.e. M K = -19.5 inside 300 pc 
radii. 

The Heart of the Bird consists of a fairly regular looking 
barred spiral (in Tf-band), with a wide opening of the spiral 
arms. The arms can be traced out to at least 2 kpc from 
the nucleus. In a radial aperture of 1 kpc the brightness is 
K=13.35 mag, ie. Mk = —23.16. In the optical, the spiral 
arms are more pronounced than the nucleus, though only 
small sections of the arms can be seen in the B-band. 

The Head, the northernmost component, appears to be 
an irregular galaxy, extending approximately 2 kpc with nu- 
merous distinctive bright knots in the K image, presumably 
bright HII regions and super star clusters (SSC). The op- 
tically brightest section is the northern part of the Head, 
while the brightest NIR knot is in the most obscured south- 
ern part of the Head. The total brightness of the Head is 
-fT=13.80 mag in a 1 kpc radial aperture {Mk = —22.71) 
i.e. only half a magnitude fainter than the Heart. Individual 
knots, excluding the more diffuse K-band background light, 
range in brightness from K = 17.6 ± 0.1 to K ~ 20.0. We 
can isolate 12 of these H II regions mainly at the edge of 
the component closest to the Heart; the blobs get fainter 
with increasing distance from the Heart and Body. Most of 
these blobs are in highly extincted regions and are not visi- 



ble in the optical images - an elaboration on point-like SSC 
candidates is presented below in Section [3. 1.41 

Additionally, there are tidal tails, the Wings, extend- 
ing East and North- West, visible for 8 kpc and 15 kpc, re- 
spectively. In addition we separate a West Wing, a shorter 
section of the tidal tail below the larger structure extend- 
ing North- West, and which appears from the ACS images 
to continue from the Heart spiral arms. All the Wings also 
have numerous bright knots though not as extended as the 
ones in the Tail. Finally, there is a very complex dusty area 
North of the E-wing and to the NE of the Head. 



3.1.2 Surface brightness profiles and effective radii 

To study the morphologies of the galaxy compo nents of the 
Bird more quantitatively, we used the G ALFIT (|Peng et al.l 
120021) software on the NACO data to fit 2D surface bright- 
ness distributions to the two brightest components, the 
Heart and the Body, and in addition fit simple surface 
brightness profiles to cuts along the semi-major axes of these 
two galaxies. All the systems of the Bird are quite disturbed, 
so unique solutions, especially with G ALFIT, are difficult to 
find. Nevertheless, results are satisfactory for the purpose of 
classifying the galaxies, and to derive effective radii to better 
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Figure 5. A Spitzer/YRS spectrum of the Bird. The main PAH features and other detected emission lines are indicated. The spectrum 
continues featureless beyond the plotted area. The small discontinuity at ~ 14 um is the transition from SL to LL mode observation. 



Table 2. Photometry of IRAS 19115-2124. Absolute magnitudes 
can be calculated with a distance modulus of 36.51. 



Component 


aperture a 


B 


I 


K 




(kpc) ~ 


(mag) 


(mag) 


(mag) 


Bird total 


15i 


15.6 


13.6 


11.0 


Body 


1.0 


19.80 


17.00 


12.82 




0.5 


22.00 


18.00 


13.58 




0.2 


> 25.7 


19.32 


15.68 


Heart 


1.0 


19.31 


16.87 


13.35 




0.5 


21.52 


18.22 


14.41 




0.2 


23.81 


19.63 


16.63 


Head 


1.5 


18.64 


16.49 


13.39 




1.0 


21.51 


17.13 


13.80 



a The 1.5, 1.0 and 0.5 kpc radius apertures are calculated above 
the global sky level, the 0.2 kpc aperture is above the local galaxy 
background. Photometric errors in all bands are dominated by 
systematic effects and are approximately at 0.05 mag level — only 
the 0.2 kpc aperture B-band values have errors of the order of 
0.2-0.4 mag. 

^ Flux within a large 30 arcsec diameter, if-band value is also 
equal to the 2MASS total magnitude 



than a factor of 2 accuracy, that is needed in mass estimates 
later. 

First of all, we do not find point-source components 
from the central parts of the galaxy system. The Heart of 
the Bird, the apparent barred spiral galaxy, is a combination 
of a de Vaucouleurs bulge component and an exponential 
disk. The best-fitting effective radius for the disk compo- 
nent is r e =1.2 kpc using a ID profile along the semi-major 
axis (see Fig. 0, and when fit out to a distance of 1 arcsec. 
GALFIT turned out to be quite sensitive to the adopted 
'sky' value, i.e. in this case the value assigned to the un- 
derlying light from the tidal tails and other neighbouring 
components. With a realistic range of experiments, GAL- 
FIT produces a range of r e — 1.1 — 2.3 kpc, consistent with 
the simpler method. Table [3] gives the best-fitting GALFIT 
parameters for all components. The derived bulge-to-disk 
ratio with these parameters is B/T — 0.09, i.e. typical of 
late type spirals. The lower left panel in Fig. [6] shows the 
residual image after subtraction of both the bulge and disk 
components using GALFIT. 

The Body of the Bird, a very disturbed galaxy, was more 
difficult to fit using GALFIT, and impossible with simple 
exponential and/or de Vaucoulers components. Using a ID 
slice through the galaxy (Fig. [7} out to a distance of 1 arcsec, 
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Table 3. GALFIT and ID profile fit parameters for IRAS 19115-2124 NACO data. 



Component 


subcomp 


Method 


K 




n 


b/a 


PA 


boxiness a 








(mag) 


(kpc) 






(deg) 




Heart.. 


bulge 


GALFIT^ 


16.3±0.2 


1.3±0.2 


4 


0.9±0.1 


-60±10 






disk 


GALFIT - 


13.7±0.2 


1.9±0.4 


1 


0.4±0.1 


7±5 






disk 






1.26±0.03 


1 




12 




Body ... 


bulge 


GALFIT 


16.1±0.2 


0.7±0.2 


2.3±0.2 


0.7±0.1 


33±10 


+0.6±0.1 




disk 


GALFIT 


12.4±0.2 


2.1±0.4 


1.6±0.1 


0.5±0.1 


4±7 


+0.1±0.1 




disk 


ID 




2.52±0.05 


2.23±0.04 




13 





a parameter c > boxy, c < disky 

^ The errors reflect systematic uncertainties of parameters found when using a range of fitting areas, initial 
values etc. Formal GALFIT errors of parameters are much smaller. 
c Formal errors are given for ID fit parameters that were not fixed. 



a single Sersic profile with n=2.2 and r e = 2.5 kpc is the 
best fit. Using GALFIT, a reasonable fit is found with two 
separate Sersic profiles with n=2.3 and n=1.6 for the nuclear 
and outer region components, respectively. Though the best 
fit here was again sensitive to the background value, and 
the index n is coupled to the radius, a reasonable range of 
parameters resulted in r e ~ 1.0 — 2.3 kpc. The right panel of 
Fig.|6]shows the original surface brightness distribution, and 
the residual after the best-fitting subtraction, which reveals 
a complicated underlying appearance, with evidence for a 
bar-like structure extending NE-SW, and perhaps remains 
of spiral arms. It is also noteworthy that regardless of the 
other parameters, the 2D GALFIT profile for the nuclear 
region of the Body is always "boxy" (with the GALFIT 
boxyness parameter c w +0.6). 



3.1.3 Search for supernovae in NACO images 

We have NACO i\-band images obtained at two different 
epochs separated by 153 days. This allows us to search for 
supernovae that have exploded within the nuclear regions of 
IRAS 19115-2124 during this period or before the first epoch 
image. The images were aligned, subtracted a nd analysed in 
the manner described in lMattila et al.l ()2007l ) . where we de- 
tected the SN 2004ip in a similar NACO dataset of the LIRG 
IRAS 18293-3413. No obvious point-sources were found from 
the subtracted images (Fig. |8j . We then estimated SN de- 
tection limits by placing simulated sources made from a real 
PSF of a nearby bright star in a grid of locations within the 
nuclear regions of the Bird prior to the image subtraction. 
We used an aperture radius of 0.15 arcsec and a sky annulus 
between 0.20 and 0.30 arcsec for measuring both the noise 
and the simulated sources in the subtracted images. This 
yielded K = 21.5 as the 3<r limiting magnitude for the nu- 
clear SN detection within the Head. Within the Body, where 
larger image subtraction residuals were seen, we obtain a 3cr 
limiting magnitude of K = 21.3. However, within the inner- 
most ~0.3 arcsec (or ~300 pc) diameter of the Body the SN 
detection is more difficult and a SN of K — 21 mag would 
remain undetected in our NACO images. 



3.1.4 Super star cluster candidates 

To create a list of SSC candidates we searched for all point- 
sources in the ACS images within the NACO image FOV 
- 83 were found - and aperture photometry was performed 
in 0.2 arcsec radii apertures and the background measured 
in the surrounding 0.1 arcsec sky ann uli. Aperture correc - 
tions for each ACS image were applied (|Sirianni et al.ll2005l ). 
and Galactic extinction was corrected for. Identical aper- 
tures were then used on the NACO image, and 78 of the 
83 point-sources were detected, 2 of them marginally (cross 
points in Figure |SJ. Aperture corrections of 0.84 mag were 
applied, determined from brighter isolated point-sources in 
the field; scatter is of the order 0.1 mag, and a small de- 
pendance on the distance to the WFS reference star was 
ignored, since it is smaller than the scatter. 

Figure[9]shows plots of the SSC candidates, the different 
symbols differentiate between point-sources inside the Bird 
system, (45 of them, open diamonds) , and those outside the 
system (31, filled circles). The largest single concentration 
of these point-sources is along the W and NW Wings. The 
number density of the measured sources inside and outside 
of the galaxy was compared to the number density of point 
sources (foreground stars) in the ACS images outside the 
NACO FOV. From this we estimate that ~ 10 per cent of 
the SSC candidates inside the Bird and ~ 90 per cent of 
those outside the system are foreground stars. Therefore, 
the colour-colour diagram (Fig. [9] lower panel) appears to 
nicely separate the SSCs from the foregound stars. A handful 
of points have very red B — I ~ 3 and I — K ~ 4 colours; 
these are all within the core regions of the Head and Heart. 
We note that based on Fig. [9] the luminosity function of 
the SSC candidates appears to be brighter by at least a 
magnitude than e.g. i n the local luminou s interacting galaxy 
prototype Antennae (I Anders et alj|2007h . However, we defer 
more detailed discussion on SSC candidates in the Bird to 
a separate work. 



3.1.5 MIR emission 

Figure [4] shows the Bird i\-band data overlaid with Spitzer 
IRAC and 24 (am MIPS contours. The 3.6 [im data closely 
matches the K-band, but the longer the wavelength, the 
more of the relative flux comes from the northernmost Head 
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Figure 6. Close-ups of Heart (left) and Body (right) galaxy nuclei of the Bird from the i\-band AO data. Top panels show the data, and 
the bottom panels shows the respective galaxy nuclei after subtracting the GALFIT best-fitting bulge and disk components - see text and 
Table[3]for the values. The brightest residuals are of the order of ~ 7 per cent of peak fluxes for both nuclei. There is a bar-like structure 
revealed in the Body nuclear region, which was not obvious in the unsubtracted image. Images are shown with inverted brightess scale. 



component, the irregular galaxy with numerous HII/SSC 
regions. We also used the 3.6 [im band map to extrapolate an 
approximate stellar contribution in the other IRAC bands 
by simply assum ing a Vega-like stellar spectrum (see e.g. 
IPahre et afll2004h . The stellar light subtracted 5.8 and 8.0 
|xm maps do not differ appreciably from the non-subtracted 
ones, which is not surprising since these are expected to be 
dominated by non-stellar emission everywhere. The stellar- 
subtracted 4.5 ixm image, however, clearly reveals the Head 
to be the major source of non-stellar MIR flux. 



MIR colours are examined in more detail in Figure [TOl 
where all IRAC pixels in a 22 by 22 arcmin area around the 
Bird are plotted. Pixels at the locations of two field stars are 
located close to zero-colours as expected. Most colours of the 
interacting system are close to those expected from late type 
spiral an d irregular galaxies, e.g. [3.6l-[4.5l~0.1 and [4.5]- 
[5. 8] ~ 1.5 (|Davoodi et alj|2006l ; ISmith et alll2007h . There are 
two significant exceptions, however. All colours with the 8.0 
(j.m IRAC band are very red, indicating stronger than nor- 
mal spiral galaxy PAH emission contribution in that band. 
Secondly, the Head (red squares) has very red colours in 
[3.6]-[4.5] and [3.6]-[5.8], more typical of e.g. Seyfert nuclei. 
This must be a result of warm to hot dust in the irregu- 
lar galaxy, which can in principle be heated either by very 
strong star formation, or AGN activity. We investigate spec- 
tral evidence for AGN further in Section \3. 2. II 



3.2 Analysis of spectra 

3.2.1 Physical characteristics from emission lines 

Several warm gas emission lines are readily detected in the 
SALT/RSS spectra, as seen for example in Fig. [2] Table [4] 
summarizes the optical line ratios available from our data for 
the different parts of the Bird, gathered from all three slit 
positions, and calculated using multi-component fitting with 
the IRAF task splot. The line ratios are not corrected for 
stellar absorption; with typical EW(ifa) > 50 A (see below) 
we expect the correction to be not significant. All areas have 
uniform ratios, in the range [N II]A6583/Ha ~ 0.40 - 0.45 
and [S H]AA6716,6731/Ha ~ 0.4 bot h suggesting HII type 
regions l|Veilleux fc Osterbrock|[l987f ). The only notable ex- 
ceptions are the nuclear region of the Heart, and the de- 
blended blue wing from the Heart: while the nucleus has even 
more solid HII type ratios than the average, the blueshifted 
emission from the Heart has an elevated [Nil] level with the 
ratio to Ha at ~ 0.6 (similar to other deblended blue-shifted 
components as well, which are not tabulated). This is close 
to the HII vs. AGN separation line; but instead of pho- 
toionization by an AGN, it may well indicate shock-heating 
mechanisms contaminating the normal ratios expected from 
HII regions. 

The strongest equivalent width of Ha in the main com- 
ponents is detected in the Tail, where EW(Ha) « 100 A. 
The next highest are in the Head and Heart. The whole mea- 
sured range is typical of HII type (U)LIRGs (|Veilleux et al.l 
Il995l) . The strongest individual star formation knots falling 
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Figure 3. The reduced 2D spectrum around the Hct and N II lines 
from three different position angles. The main components of the 
merging system are indicated, see Fig.fl] The bright emission lines 
can be traced for over 20 arcsec, i.e. 20 kpc. 



within our slits also result in EW(Ha) ~ 100 A values (such 
as a blue double-knot in the W-wing 8 arcsec West of the 
Heart nucleus) . The EW values imply ages of 6 to 7 Myr for 
stellar populations involved, with an instantaneous bur st, 
and depending on metallicity (|Monreal-Ibero et al.ll2007h ■ 

A rough averaged metallicity estimate can be ob- 
tained usi ng the N2 = [N II] A6584/Hq r atio: from re- 
lations in iDenicolo. Terlevich fc Terlevichl (|2002T ) we get 
12+log(0/H)= 8.85 ± 0.12 where uncertainties both in the 
ratio and the method are included. The metallicity would be 
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Figure 7. i^-band surface brightness profiles of the smaller 
barred spiral (Heart, left), and the larger NIR galaxy (Body, 
right). The dots represent all pixel values as a function of dis- 
tance from the nuclei, and red squares are values averaged over 
0.1 arcsec bins along the major axis only. The dots above the 
squares in the left panel come from the bars. The blue lines are 
chi-square fitted profiles along the major axes: an exponential and 
de Vaucouleurs component in the left panel, and a single Sersic 
profile on the right, see Table [3] for best-fitting values. 



Figure 8. The two epochs of NACO K-b&nd data subtracted 
from each other (Sep 10 - Apr 13). No point-sources, SNe, are 
found down to K = 21.3 mag. The image shows mostly ran- 
dom noise but some significant residuals, due to differences in the 
seeing halos between the two images, are apparent around the 
location of the bright Body of the galaxy. The brightness scale is 
inverted and the image area is identical to those in Fig. [4] 
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Figure 9. SSC candidates in the NACO image FOV, i.e. all point- 
sources within the galaxy system (diamonds) and outside it (cir- 
cles). Triangles show detections within core regions of the Head 
and Heart and crosses denote marginal NIR detections. A red- 
dening vector of E(B — V) = 0.32 is indicated. 



8.75 for the lowest detected ratio at the deblended nuclear 
component of the Heart. 

The R [SII ] = [SII]6716/[SII]6731 ratio ranges between 
R[ SII ] « 1.0 — 1.5 throughout the Bird, which correspond to 
electron densities from n e = 500 cm -3 to less than 10 cm" 3 . 
Apertures within the Heart and Body show the highest den- 
sities of 300 — 500 cm -3 , while all the other Bird components 
appear to have n e < 20 cm -3 . 

From the IRS mid-IR spectrum (Fig. [5} we measure a 
line ratio of [NeIII]15.5n/[Ne II] 12.8^=0.10, which is typical 
of low to medium excitation starbursts, and corresponds to 
effective stellar temperatur es of 3.5 — 3.7 x 10 K in H II 
regions (|Giveon et alj|2002h . The IRS spectrum is examined 
in more detail in Section f4. 31 



3.2.2 Velocities of major components 

As seen from the RSS 2D spectra (Fig. [3} the velocity 
structure of the Bird is quite complex, and we used multi- 
component fitting to derive line ratios (above) and accurate 
velocities and dynamics (next and later sections). However, 
we found it useful to also simply fit a single-component gaus- 
sian to the Ha lines along all the three slits for general char- 
acteristics of the primary components. The resulting velocity 
curves, each point representing a single pixel row, are shown 
in the bottom panels of Fig. 1111 where the main structural 
components are also indicated. The results are unchanged, 
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Figure 10. Pixel-by-pixel IRAC colours of IRAS 19115-2124 
area. Black crosses are pixels within the the galaxies (defined as 
/8.o > 6 MJy/sr), and the green triangles are the two bright stars 
in the field. The blue diamonds pick out the area of the NW wing, 
and the red squares the Head of the Bird - the latter is clearly 
separated with the [3.6]-[4.5] colour. The tidal tails, in contrast, 
do not differ from the bulk of the system. 



though more inaccurate, if the weaker [Nil], [S II], or [OI] 
lines are used. 

The Body and the Tail of the system join together 
smoothly in velocity space, while the Head is clearly a sepa- 
rate component from the rest. The Heart, the smaller spiral, 
is the only component that shows a clear rotational structure 
with peak-to-peak velocity of approximately 370 km s _1 , 
while any rotational velocity of the other components is 
< 100 km s _1 from our slit orientations. The velocity of 
the W Wing appears to join well the arms of the small spi- 
ral, but there is a jump of 100 km s _1 from the arms on the 
other side to the E wing. 

A general picture emerges from these spectral profiles, 
where the Head is moving away from us with the highest 
radial velocity of 14 940 km s _1 , i.e. +360 km s" 1 with 
respect to the systemic velocity (see below), and the E Wing 
has the second highest relative velocity of +240 km s - . 
The NW Wing, in contrast, is turning our way with a high 
relative velocity of —400 km s~ . The Tail has an average 
of —110 km s _1 offset with respect to the systemic velocity, 
while the Body appears to be very close to the systemic 
velocity. 

Velocities of th e different component s have an excellent 
match with those of lMirabel et all l|l990l) . whose CO(l - 0) 
based values show three distinct components at 14400-14550 
km s" 1 , 14600-14800 km s"\ and 14850-15000 km s~\ cor- 
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Table 4. Properties of the components of IRAS 19115-2124. 



Component 


subcomp 


R a 
(k~) 


(km s — 1 ) 


D^~+ ( lo^}^ 
u VOt \ 1^— 

(km s — 1 ) 


(km s -1 ) 


[NII]/Ha e 


[SII]/Ha/ 


EW(Ha|£ 

(A) 


Heart... 


totaljj. 


4 


14520±20 


185±6^ 


129±7 


0.43 


0.35 


51 




nuclear-? 


0.5 


14576±9^ 




76±9 


0.32 


0.28 


35 




blue wing^ 




14362±23 






0.56 


0.39 


19 




W-disk^ 


1.5 


14362±12 




90±6 


0.44 


0.41 


37 




E-disk^ 


1.5 


14701±6 




52±5 


0.42 


0.42 


39 


Body 




3 


14590±40 


75±4 


134±16 


0.42 


0.32 


42 


Head 




2.5 


14940±30 


40±5 


93±8 


0.47 


0.45 


69 


Tail 




6 


14467±7 






0.41 


0.40 


97 


E-wing 




8 


14817±6 






0.45 


0.40 


28 


NW-wing 




16 


14175±6 






0.44 


0.45 


27 


W-wing 




3 


14337±6 






0.38 


0.34 


63 



a Approximate extent of component in the NACO image, not formal r e 

b Heliocentric line-of-sight velocity; error includes formal error and the spread of velocities derived from different lines 

c Half of peak-to-peak velocity from rotational curves of one-component fits, not corrected for inclination effects 

" Velocity dispersion, corrected for instrumental resolution and redshift 

e From [NIIJA6583, not corrected for reddening, uncertainties approximately ±0.03 

f From [S II] AA6716+6731, not corrected for reddening, uncertainties approximately ±0.03 

3 Uncertainties approximately ±3 A 

^ Wide aperture, unblended multiple components 

1 The deblended adopted rotational velocity from this aperture is 157 km s — 1 . See text. 

■? Deblended component 

"* Adopted systemic velocity of the Bird 

' Deblended blue component of the nuclear aperture 

m Separate aperture West of the nucleus 

n Separate aperture East of the nucleus 



responding to our Tail, Heart, and Head, respectively. More- 
over, their fourth, weaker and narrow component at 14200 
km s -1 agrees exactly with our derived velocity for the NW 
wing of te Bird (see Fig. I12|l . 

We also plot the velocity dispersions (middle panels 
of Fig llip extracted from the one-component fits along the 
slits. The Wings and the Head are kinematically cold, while 
the Heart and Body nuclei have quite high values (we return 
to velocity dispersions for the purpose of mass estimates in 
Section l"3.4p . Note also that the strength of the Ha line emis- 
sion (top panels of Fig. Hip does not always exactly corre- 
spond to the position of the maximum velocity dispersion. 
The velocity dispersion peak does, however, match well the 
position of the peak continuum emission along the slit, which 
is chosen as the zero-point of the x-axis radial distance. 

Table U summarizes the properties of the major compo- 
nents of the Bird. 



3.2.3 Velocities of sub-components and the systemic 
velocity 

The lower panel of Figure [13] shows the spectral profile of 
the October 22 data, at the wavelength region of Ha and 
[Nil] lines, over a one-pixel (0.26 arcsec) aperture at the 
centre of the continuum of the observed spectrum, which we 
expect to correspond closely to the centre of the nucleus of 
the galaxy in the NIR image. A two-component gaussian is 
fit to the Ha profile, the heliocentric velocities of which are 
plotted in the figure. Exactly the same velocity components 



were fitted for the [NII]A6583 line. The plotted [NII]A6548 
curves, in contrast, are not fits, but are merely one third of 
the strengths of fits to the [NHJA6583 line, along with the 
appropriate wavelength shift. 

The upper panel of Figure [13] shows the ID spectrum 
profile of the same slit and same location, but summed in- 
side a wider aperture of 1.5 arcsec. Wider apertures give 
essentially the same result. Comparing to the NACO image 
(Fig. [1} it is clear that the aperture should include a signif- 
icant part of disk and spiral arms, in addition to the cen- 
tral bulge. The main features of the Ha are now fitted well 
with three gaussian components. The central component is 
kept fixed at the same velocity and width as derived from 
the small aperture profile. A higher velocity component not 
present in the small aperture profile is evident, and the lower 
velocity component has become stronger. The same velocity 
components were again fitted for the [NII]A6583 line, and 
[NII]A6548 was treated as before. 

The wider aperture profile shows a typical double-horn 
structure of spiral galaxies, with an additional blended cen- 
tral component. This latter middle velocity component of 
14576 ± 9 km s _1 must be from the nuclear region, since 
it is present identically also in small apertures. The value 
is close to the average of the low and high velocity compo- 
nents, i.e. the two sides of the disk. Independent of these 
spectral fits, we also derive a value of approximately 14560 
km s^ 1 for the point of symmetry in the rotation curve. 
Since the Heart of the Bird is the central component of the 
whole system, and since this velocity corresponds spatially 
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to the AT-band continuum peak, we adopt this 14576 ± 9 
km s -1 as the systemic velocity of IRAS 19115-2124. 

Note that the strongest Ha flux does not coincide spa- 
tially with the location of systemic velocity, i.e. the contin- 
uum peak position (see right-most panels of Fig. Hip . We 
measure v = 14622 ± 9 km s" 1 for the deblended velocity 
at the strongest Ha location, approximately 1.0 arcsec East 
of the nucleus. Without the high-resolution AO image, this 
would have naturally been interpreted as the central veloc- 
ity. Indeed, the literature systemic velocity of 14608 ± 48 
km s" 1 (NED, IStrauss et all Il992r ) agrees well with this 
strongest Ha component. It should be noted, however, that 
none of these values agree with the NED values for individ- 
ual IR AS 19115-2124 N and S components from iKim et all 
1 19951 ). 

Returning to Fig. 1131 upper panel, there are some ad- 
ditional velocity components we did not fit in the Heart 
central regions. In the wide aperture extraction there are 
two weaker structures at 14820 and 14950 km s _1 which are 
contamination from the Head and the E-Wing. But more 
significantly, the blue side of the lowest velocity (vl) com- 
ponent is not well fitted: the FWHM=3.lA of this compo- 
nent is wider than the 2.7A of the higher velocity arm, al- 
ready implying contamination, and the fit still leaves excess 



emission. Moreover, the blue excess wing is present also in 
the small aperture profile (lower panel), i.e. in the nucleus 
only aperture. We interpret these effects as outflowing gas 
present in most apertures over the system, blue-shifted here 
by velocities exceeding 150 km s~\ That the blue compo- 
nent has an elevated [N II] A6383/H q ratio may be a result 
of contribution from shock -heating (|Veilleux fc Rupkell2002l ; 
iMonreal-Ibero et aLlliooih . 



3.2-4 Outflows detected in Na D absorption 

Additional evidence for the outflows come from strong Nal 
D AA5890, 5896 (NaD) absorption doublets in the spec- 
tra, the strongest case of which is seen in the central 
aperture of the Heart shown in Fig. 1141 In strongly star- 
forming systems, the NaD absorption originates partially, 
or mostly, from ISM, rather than being only of stellar ori- 
gin, and it has often been used specifically to st udy out- 
flows from star-forming galaxies and (U ) LIRGs ([Phillips! 
| 1993 |: ISparks. Carollo fc Macchettol 1 19971: iHeckman et al.l 
20od : [Rupke, Veilleux fc Sanders! |2002| ; ISchwartz fc Martini 
booi Martini I2005T I. and to isolate kinematics of different 
gaseous and stella r components with in the galaxies (e.g. 
lArribas et alll2007l ; iBrosch et al.ll2007l ). 
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Figure 12. We have extracted the CO(l-O) observation of 
IRAS 19115-2124 from of Mirabel et al. (1990) and labelled in 
the figure the velocity components which match the optical spec- 
troscopy derived velocities from this work. Note also the similar 
strenghts of CO emission, and hence molecular hydrogen mass, 
of all three components. 
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Figure 13. Extracted spectral profile (black solid line) at the Ha 
and N II region of the PA=90 spectrum. The upper panel shows 
the profile over a wide 1.5 arcsec aperture at the center of the 
Heart galaxy, while the lower panel is a single pixel, 0.26 arcsec, 
aperture at the same location. Blue dotted line shows the individ- 
ual fitted components, the red line the total fit. The [NIIJA6548 
profile is affected by a telluric absorption band at 6870A. 
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Figure 14. A three-component fit (red solid line is the total of 
dotted blue individual fits) to the Nal D AA5890,5896 absorption 
doublet and the He I emission line in the PA=90 slit at the nuclear 
region of the Heart. The derived velocity of the absorbing com- 
ponent, 14400 km s — 1 , fits well the excess blueshifted emission in 
Ha and [Nil] lines (Fig. 1151. 



The measured velocity of 14400 km s _1 for the doublet 
coincides with the blue wing of the nuclear aperture profile 
of the Ha and [N II] lines, implying an outflow velocity of 
the gas of approximately 180 km s _1 at this location. 

Most slit positions and apertures show NaD absorp- 
tion, usually blueshifted from the emission line positions. 
The doublet profile is not always as clear as in the case 
of the Heart, but rather shows a wide range of velocities, 
and/or the presence of a wider absorption component pre- 
sumably originating from zero-velocity stellar component. 
Table [S] summarizes the outflow characteristics in the var- 
ious components of the Bird, based on NaD. Where possi- 
ble, we have also extracted the equivalent width ratio Rn & 
= EW (5890) / EW(5896), which is a measure of optical 
depth (|Rupke. Veilleux fc Sandersl20o3 ; ISchwartz fc Martini 
|2004| ): R,Na = 2 for optically thin medium, while the ratio 
approaches unity with infinite optical depth. 

Overall, the strongest outflows are related to the loca- 
tions of the Heart and Head, and they are somewhat weaker 
over the Body. If the outflows were from the major galaxy 
nuclei, the Heart and Body, the implied radial velocities are 
in the range 100-300 km s~ . However, the velocities can be 
up to 600 km s _1 if some of the outflowing material origi- 
nates from the Head. Any outflows at the positions of the 
Wings or the Tail are below our detection limits. 



3.3 Extinction estimate and gas and dust content 

Since our RSS spectra did not reach the H/3 line, we first esti- 
mate an average extincti on from the emiss ion line strengths 
of IRAS 19115-2124 from lKim et"aH l| 19951 ): their Hcv/H/3 » 
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Table 5. Outflows and related velocities from IRAS 19115-2124 
components. 



Component NaD doublet^ NaD range_ Optical depth^ 
(km s _1 ) (km s —1 ) 



Heart 180 180 - 40o£ thick 

Body - 270 thick 

Head 300 - 600 
Tail 

E-Wing thin 

NW-Wing 

W-Wing 



and fcont are the observed flux density and the estimated 
continuum flux density, 5*sil ~ —1-4, implies optical depths 
in the range ry = 30 — 200 depending on the dust shell 
thickness and geometry. This range can be constrained 
by making use of the rest wavelength continuum ratio, 
/„(14p.m)// l ,(30nm) ss 0.15; the geometric models giving 
consistent optical depths with both indicators simultane- 
ously are ones with very thick Rout/Rin > 200 and nearly 
uniform radial density distribution (p ~ for r _p ), which 
indicate extinction at the lower end of the range, tv ~ 30 
(Ay ~ 30). Thus, the MIR spectrum extinction indicators 
show the Bird to be a moderately extincted source as far as 
(U)LIRG nuclei are concerned. 

Finally, we make an estimate of the gas masses involved. 
Taking an optical extinction estimate of Ay = 4 and using 
a standard gas-to-du st relation of N(H)/Ay ~ 1.9 x 10 21 
atoms cm -2 mag -1 (jBohlin. Savage fc Drake] 1 19781 ) we get 
n(H) = 1.3 x 10 -5 kg cm -2 , so assuming an absorption layer 
covering the whole system, say, a 3 by 8 kpc area, the gas 
mass would be M ~ 1.5 x 10 9 M@. If, instead, the higher op- 
tical depth of Ay = 30 were used in conjunction with a 2 kpc 
diameter area of the 8 and 24 [im emission around the Head 
and Heart components, we would arrive at M ~ 2x 10 9 M . 
These are reasonable values rec alling that the CO m easure- 
ment of molecular hydrogen bv lMirabel et all (|l990h . prob- 
ing the full emitting volume regardless of geometry and op- 
tical depth, gave 3 x 10 10 M . 



3.4 Masses and dynamics 

3.4-1 NIR light based mass 

Before estimating the dynamical mass of IRAS 19115-2124 
and its components, we use the absolute K-band magnitudes 
presented above in Section 13.1.11 The total integrated ab- 
solute magnitude from our NACO image is Mk = —25.5 
mag, or about 3.6 x L*. The individual brightnesses in 
1 kpc aperture radii for the Body, Heart, and Head are 
M K = -23.7, -23.2, and -22.7 mag, respectively. These 
translate to 6.3, 4.0 and 2.5 x 10 10 Lk,q (or 0.7, 0.4, and 
0.3 L* , and nearly 2/3 of the total Bird brightness is thus 
contributed from outside the 1 kpc apertures). It is com- 
monly assumed that NIR light traces reasonably well stel- 
lar mass (though the quantitative relation is wrought with 
complications) - as an estimate, we use a if-ba nd mass-to- 
light ratio from iThronson fc Greenhouse! (Il98cf ). and calcu- 
late 3.0 x 1O U M for the whole Bird system (if =11.0 mag, 
and D L =207 Mpc). Taking m* = 1.4 x lO n M , this corre- 
sponds to a mass of 2.1 m* in (old) stellar population in the 
Bird. The three individual components, within 1 kpc radii 
each, would range from 0.16 m* (Head) to 0.40 (Body), 
see Table [6] Since the component masses were calculated 
within 1 kpc radii the likely total masses would perhaps be 
a factor 2 larger, in accordance with the total-light derived 
value. 

No extinction corrections were made above, making 
these estimates lower limits. Average extinctions over the 
major nuclei are in the range Ay ~ 2.5 — 5, as seen above, 
which would mean upward mass corrections of 30 to 60 per- 
cent. However, we have not attempted a correction due to 
light from bright, recently-formed red supergiant popula- 
tions, the continuum light of which would in contrast result 



a Offset velocity of the NaD doublet relative to the component's 
velocity from emission lines. '0' means a zero-velocity feature, 
and '-' a non-detection. 

^ Velocity range of a wider NaD absorption feature. 

c From the R Na = EW(5890)/EW(5896) ratio. '-' means the 

ratio can not be measured reliably. 

" In a wider aperture than shown in Fig. 1141 

9.C0, implying Ay « 3.3, or Ak ~ 0.36, which is typical of 
local ULIRGs. Of course, it is obvious from the optical im- 
ages that extinction is highly variable over the Bird, for ex- 
ample the Body and southern parts of the Head are most af- 
fected by dust, but as an averaged value the above indicates 
significant, yet not extreme extinction. This value is con- 
sistent with a statistical value found from an empirical rela- 
tion fou nd for IR-luminou s galaxies in the Spitzer First Look 
Survey (|Choi et al J 1 20061 ): Ay = 0.75 x \og(L IR /L @ ) - 6.35 
which would give Ay = 2.5 for the Bird. 

We also make an extinction estimate based on the de- 
tected NaD absorption feature - this can be done assum- 
ing it is most l y of i nterstellar origin, rather than stellar. 
IVeilleux et aT] jl995h show that there is a correlation be- 
tween EW(NalD) an d colour excess E(B-V) for HII typ e 
(U)LIRGs (see also ISparks. Carollo fc Macchettol Il997l) . 
Taking EW(NalD) ~ 5x E(B-V) and our EW(NaID)« 4 A 
from the strongest case in the Heart, gives E(B-V) ~ 0.8, 
and hence Ay ~ 2.4. 

Broad-band colours give indications of dust redden- 
ing in specific regions. Using photometry in 500 pc radii 
for the Body and Heart nuclei (Table [5]) and in 1 kpc ra- 
dius for the brightest NIR feature in the Head, we measure 
I — K = 4.4, 3.8, and 3.7, respectively, and B - K = 8.4, 
7.1, and 6.6. Since normal spiral galaxy colours are approx- 
imately B — K ~3.5 and I — K ~2, these colours imply 
Ay ~4-5, i.e. somewhat higher extinctions than from the 
methods above. Both the Body and Heart regions would 
in fact nominally be classified as E xtremely Red Objects, 
EROs i|Vaisanen fc Johanssorj|2004h . 

The Spitzer d ata has the poten t ial of probing the most 
extincted regions. iLevenson et al.l (|2007l ) have calculated 
models that relate optical depths to IRS spectral features: 
the strength of the MIR silicate feature (see Fig. [SJ) de- 
fined as Sstl = In [/ o6s (9.7|i.m)// corl t(9.7|J.m)], where f ob3 

3 The value is from their North component of the galaxy, since 
their South component clearly suffers from a simple line misiden- 
tification. 



16 P. Vdisdnen et al. 



in us predicting too much mass: their contributions range 
from nearly negligible in normal spiral galaxies to perhaps 
40 per cent in localised youn g starburst regions fe.g. lRhoadsl 
Il998l ; [james fc Seigar][l999h . The correction may be signif- 
icant in case of the Head, where intense individual (pre- 
sumably young) clusters are observed in the NIR, but it is 
unlikely to be very significant in the Heart and Body com- 
ponents. As a further pot ential correction, ba sed on typical 
NIR spectra of ULIRGs (|Murphv et al.lll999h . we estimate 
the expected contribution from emission lines (especially 
Paa, Br7, and H2 vibrational series) to the broad-band light 
to be maximally ~ 10 per cent, but likely much less than 5 
per cent, since the strongest contamination from Paa only 
marginally enters the edge of the NACO K 3 filter at the Bird 
redshift. Taken together, it is more likely that an upward 
correction to mass due to extinction dominates over any 
downward corrections due to starbursts and emission lines 
from gas. Since all the corrections are somewhat uncertain, 
but nevertheless of similar magnitude, we prefer to keep the 
_K"-band light mass estimate as is, as a lower limit, with- 
out applying the corrections. Finally, note that the mass- 
to-light conversion used implies values of M/L~ 0.9 in the 
system: this is consistent with values expected from stellar 
populations in high-luminosity (isolated) spirals calculated 
from models that also include moderate starburst episodes 
jBell fc De Jondl200ll ). 

3.4-2 Dynamical mass 

For dy namical mass estimates w e follow iDasvra et al.l 
l|2006bh (see also lColina et~ai1l2005h and express m (in units 
of Mq) as 

m = 4.7x 10 5 (3cr 2 +tv ot ) r e (1) 

where a and v ro t are the line-of-sight velocity dispersion and 
the inclination corrected rotational velocity, both in units of 
km s _1 , and r e is the effective radius in kpc. The relation 
assumes a virialised central region, and the constant term 
depends on the distribution of matter in the system: the 
value above assumes a King model with constant mass-to- 
light ratio, with tidal-to-core radius ratio of 50, which is 
midway between values for dwarf and giant ellipticals. Note 
also that the rotational velocity component effectively adds 
the mass of any disk to that of the bulge component esti- 
mated by the velocity dispersion term. 

To get the true rotational velocity v rot one has to cor- 
rect it both for deviations of the slit orientation from the 
major rotational axis and inclination effects. In the follow- 
ing v ro t(los) is the line-of-sight rotational velocity from our 
data, 9 is the angle between the slit and the major axis of 
rotation and i the inclination. 

Vrot = v TO t{los) /(cos(0) sin(i)). 

The only case of significant rotation seen in the three slit 
positions, is the Heart. The best-quality slit position PA=90 
accross the Heart is aligned nearly perpendicular against the 
major rotation axis, in principle making the required orien- 
tation correction difficult. However, the 1.5 arcsec slit-width 
actually encompasses a large fraction of the disk area, which 
we estimate from the NACO image as approximately 2.0 arc- 
sec across, and therefore we can estimate which fraction of 
the maximum rotational velocity enters the slit. Based on 



this geometrical assessment, we adopt an approximate value 
of 9 as 30 in the above equation. We adopt v ro t(los) — 157 
km s _1 measured from the fitted blue and red velocity com- 
ponents in the central 1.5 by 1.5 kpc region of the Heart 
(Fig. [13] We thus have v rot (obs) = 181 km s _1 as the ob- 
served rotational velocity. For an inclination estimate we use 
the elongation of the Heart disk, which has the ratio of the 
minor-to-major axes b/a m 0.75 from NACO isophotes, i.e. 
i ~ 44deg, assuming 0.3 for the ratio of thickness and trun- 
cation radius of the disk. Thus, v rot = 261 km s" 1 , which 
is a reasonable value for fairly massive disks. The rotational 
velocity would be v ro t = 226 km s _1 if only the inclination 
correction is applied. 

Similarly, we estimate v rot = 75 km s _1 for the Body, 
and v ro t = 40 km s _1 for the Head, both from the PA=192 
data. For the Body 9 m 0, and inclination i w 60 were used, 
whereas no corrections to the apparent observed v ro t{los) of 
the Head were made due to irregularity and ambiguity of its 
shape. 

We next measure velocity dispersions in 1.0-1.5 kpc 
apertures at the centers of the galaxies, corresponding to 
their adopted systemic velocities, with multiple compo- 
nent Gaussian fits to Ha, [NIIJA6583 and [S II] lines, and 
quadratically correct them for our instrumental resolution, 
and adjust for redshift (Table [S]). Single component fit re- 
sults, which may be contaminated by gas outflows, are also 
tabulated. The Heart and Body velocity dispersions come 
out to be o = 76 i 9 km s _1 and a = 94 ± 10 km s _1 , re- 
spectively, which can be considered as definite lower limits 
for velocity dispersion. 

Now, with the values of rotational velocities and veloc- 
ity dispersions at hand, we may estimate the masses using 
Eq.[T](see Table H)}. For the sizes r e we trust the fits over ID 
cuts along major-axes more than the GALFIT best-fitting 
values, which are very sensitive to e.g. choices of sky back- 
ground and exact Sersic n parameter. Half-light radii mea- 
sured photometrically from the NACO image are typically 
close to the lower range of r e values of both the ID cuts and 
GALFIT for the Body and Heart. We conservatively esti- 
mate the uncertainties to be at about 30 per cent level. In 
the case of the Head, the irregular morphology prohibits a 
formal r c fit, and the value tabulated is an estimate based 
on half-light radius from photometry. 

The dynamical mass of the Heart, using Eq. [T] comes 
out to be 0.36 m* (0.50 to*) using the deblended (unblended) 
velocity dispersion. Since this is a clear case of a spiral 
galaxy, we also tabulate the mass obtained directly from 
a simple Newtonian system enclosing a mass inside 3 kpc 
radius and rotating at 261 km s _1 . This estimate gives a 
lower limit of 0.32 m* for the Heart system, with the caveat 
that some part of the adopted v ro t may be tidal rather than 
orbital motion. The dynamical mass estimate of the Body is 
0.27 m* (0.50 TO*). This system is dynamically hot - ignoring 
the rotational component does not appreciably change the 
mass estimate. Finally, the corresponding dynamical mass 
for the Head is 0.08 to*. 

The mass ratio of the three main components of the 
Bird thus becomes 4.5:3:1 (6:6:1) for Heart:Body:Head when 
the deblended (unblended) velocity dispersions are used. If 
the rotational velocity estimate only is used for the Heart, 
the ratio is 4:3:1, and when if-band light based estimate is 
used, the ratio is 1.5:2.5:1. This last ratio is the only one 
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measured within identical radii. Every option nevertheless 
implies two close to equal major merger components, and a 
third component in the range 20 to 40 per cent in mass of 
the more massive ones. 



4 DISCUSSION 

4.1 Dynamics and kinematics 

4- 1.1 Reliability 

Can masses determined from velocity dispersions and ro- 
tational velocities calculated from optical emission lines be 
trusted for dynamical mass estimates? The answer depends 
very much on how closely the velocities and structures of 
warm ionized gas are coupled to the stellar stuctures and 
velocities. From the outset, there is no reason to expect 
closely matching correlations since the optical emission lines 
trace regions well outside the more extincted nuclei and since 
the very clumpy and disturbed nature of nuclear regions in 
(U)LIRGs m ake correspondences highly complicated. 

Recentlv lColina et all (|2005h presented a detailed study 
of a sample of 11 ULIRGs with angular resolutions of ap- 
proximately 1 kpc throughout, and covering a range of 
ULIRG types, utilizing optical integral field spectroscopy 
(tracing the gas distribution with Ha light), medium- 
resolution Keck and VLT near-IR spectroscopy (tracing stel- 
lar light with CO absorption bands), and millimetre CO ob- 
servations (tracing the cold gas structures). They find, as 
expected, that in general over a few kpc the ionized gas has 
a complex velocity structure, with the peak-to-peak velocity 
differences, and often also velocity dispersions, dominated by 
tidal tails and tidally induced flows, rather than rotational 
systems (though they do find ionized gas matching stellar 
velocities much more accurately than the molecular gas). 
However, they also show that the central warm gas velocity 
dispersion correlates well with the stellar velocity dispersion 
(they find a ratio of 1.01±0.13) and conclude that the central 
ionized gas velocity dispersion is, in fact, a robust tracer of 
dynamical mass of these systems, provided the emission lines 
in question can be spatially associated with the true nuclei 
identifiable virtually onl y by high-resolution near-IR imag- 
ing. We also note that iRupke. Veilleux fc Sanders! (I2002T I 
found a good correlation between the widths of emission 
lines, the Ca II triplet widths, and of the absolute K-band 
brightness of their sample of ULIRGs. 

During t his study we have s een exactly the problems 
posed in the IColina et alj (|2005h work - without the ex- 
tremely high quality near-IR imaging it would have been 
impossible to determine which features in the optical spectra 
correspond to nuclear and surrounding locations. For exam- 
ple, we are now able to tie the systemic velocity to the nu- 
cleus of a major component, rather than a bright off-centre 
HII region where the literature value seems to have origi- 
nated from. We have furthermore been very conservative in 
deblending velocity structures in spectral profiles within the 
NIR AO-determined nuclear regions to account for outflows 
and thus avoided over-estimating masses by factors of 2-4 
because of them. 

We thus conclude that the optical spectroscopic method 
of defining kinematics and dynamics of the Bird system, in 



conjunction with NIR AO imaging, should be robust and 
reflect true dynamical masses within observational errors. 

4-1.2 Tully-Fisher relation 

The Heart spiral galaxy has an absolute brightness of Mk = 
—23.2 and its rotation velocity is Vmax = 260 km s _1 . Plac- 
ing these values into a _ff-band Tully-Fisher relation diagram 
|Conselice et alj|2005h . we see that it lies well inside the 
expected 3a region of normal disk galaxies, on the fainter 
and/or faster rotation side of the relation; the small offset 
might be partially explained by under-estimated brightness 
due to extinction, or perhaps by dynamical effects due to 
the on-going interaction. 

4-1.3 The Fundamental plane 

Averaged within the respective half-light radii of the Heart 
(1.2 kpc) and Body (2.5 kpc) components, we calculate K- 
band surface brightnesses of 14.8 mag arcsec -2 and 15.8 
mag arcsec -2 . Positioning these on the Fundamental plane 
of elliptical galaxies, we find them both to occupy typi- 
cal regions of (U)LIRGs l|Genzel et al l l200ll ; iDasvra et all 
2006b), i.e. well within the scatter of the relation over dif- 
ferent types of objects, but typically 1-2 mag brighter than 
ordinary hot E/S0 galaxies. That our nuclei, together with 
other (U)LIRGs, occupy regions relatively close to the Fun- 
damental plane is commonly interpreted as evidence that 
these objects are ellipticals in formation. 

4-1-4 Escape velocities 

Following IRupke. Veilleux~fc Sanders! l|2002ri we define an es- 
cape velocity v esc at radius r for a singular isothermal sphere 
truncated at r max as 

v esc = V2v c [l + In (r max /r)] 1/2 

where v c is the rotation speed of an object in orbit. The 
escape velocity is not very sensitive to the choice of r max /r, 
though this ratio is likely to be large (> 10) since the dark 
matter halo extends far from the optically observable galaxy. 
Choosing r max /r = 10 and using v c « 260 km s" 1 (Heart 
rotation) as the extreme case, we get v esc ~ 670 km s _1 . If 
the radial velocities of the components are close to the true 
relative velocities, it would be unlikely that any significant 
fraction of the outflow material or sub-systems are escaping 
into the intergalactic medium. The 350-400 km s _1 relative 
radial velocity offsets of the Head and the NW Wing are, 
however, sufficiently close to the v esc that it is difficult to 
say whether they are truly bound or not. If one considered 
only the velocity dispersions in the two nuclei, which are at 
most « 130 km s _1 , the escape velocity would be v esc = 330 
km s _1 meaning that the Head, Wings, and much of the 
outflows would be heading out of the system based on radial 
velocities only. 

4.2 Classifying the Bird sub-systems 

4-2.1 Differentiating nuclei from star clusters 

The Bird consists of two massive components, one still a 
clear spiral, and the other a disturbed disk/bulge system, of 
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Table 6. Dynamical parameters and mass estimates for IRAS 19115-2124. We adopt m t = 1.4 X 10 1:l Mq. 



Component 


note 


M K 
(mag) 


r°_ 

(kpc) 


Vrot 

(km s -1 ) 


<7 

(km s — 1 ) 


Vrot /a 


miighl°. 

(m*) 


(m*) 


Heart 


d 








76±9 


3.4 


0.25 


0.36 




e 


-23.16±0.07 


1.26 


261±35 


129±7 


2.0 


0.25 


0.50 




I 




3.0 


261±35 






0.25 


0.32 


Body 


d 








94±10 


0.8 


0.40 


0.27 




e 


-23.69±0.07 


2.52 


75±7 


134±16 


0.6 


0.40 


0.50 


Head 




-22.71±0.07 


0.84 


40±5 


93±8 


0.4 


0.16 


0.08 



a We estimate the systematic uncertainty to be m 30 per cent in the adopted effective radii. 

k NIR light measured in 1 kpc radius, see text for M/L ratio, which also dominates uncertainty. 

Uncertainty is m 40 per cent, dominated by the uncertainty in r e 

^ Deblended velocity dispersion 

e Single component velocity dispersion 

f Using Keplerian rotation at 3 kpc 



two tidal tails/wings, and an 'extra', but rather massive, ir- 
regular component with many individual bright star-forming 
regions. 

Though there have been reports of a significant frac- 
tion of (U)LIRGS being mergers of multiple components 
|Borne et all l200fj| ; ICui et all [20011 ). most surveys tend to 
agree that they usually are interactions and mergers of two 
galaxies, wi th the incidence of multiple mergers less than 
5 per cent (Ve illeux et ail l2002l ; iBushouse et all [2002). The 
complicating factor, acknowledged in all (U)LIRG studies, is 
that any studies made in the optical regime, even at I-band, 
suffer from very complex extinction effects making photo- 
metric or morphological identification of nuclei difficult. In 
the case of the Bird, if only the /-band were considered, we 
would have 6 bright clumps (measured in 600 pc diameter 
radii) of Mi < — 17 mag, which was selected as the cut-off 
between l egiti mate n uclei and giant star-forming regions in 
ICui et al.1 (e.g. l200ll) . Only our NACO image reveals that 
three of those clumps are in fact the nucleus and two spiral 
arms of the same galaxy, and one of the clumps is a bright 
knot within the tail of the Bird, clearly a minor feature in 
the system. 



4-2.2 The nature of the Head of the Bird 

It is the more irregular feature, the Head of the Bird, which 
we are now concerned with. Is it a legitimate nucleus, a 
remnant of a pre-existing galaxy, or merely a collection of 
star clusters, formed out of the material of the other two 
main components during the interaction, for example a tidal 
dwarf galaxy (TDG)? 

It is clear that the differences in brightness are not that 
large between the three components. In the B-band and 
I-band, the Head is in fact the brightest single feature in 
the whole interacting system. Even in the if-band, which 
more closely follows the underlying mass, the Head is only 
1 mag fainter than the brightest component of the whole 
system. The dynamically derived mass for the Head ranges 
from one-third to one-sixth of the most massive component, 
de pending on the meth od used. The CO-line measurements 
of lMirabel et~ai1 (|l990h indicate an H2 mass for the Head 



of not less than half of the most massive H2 component 
(Fig. ll2[) . Finally, even the brighest NIR knot alone, which is 
totally obscured in the .B-band image, has Mk ~ —20.4 and 
Mi ~ —16.7, making it very close to the e.g. the lCui et al.l 
|200ll ) criterion for individual nuclei. 

Furthermore, the Head is kinematically totally separate 
from the Heart and Body structures, with a difference in ve- 
locity nearly 400 km s _1 (see Fig. Illf) . It also has evidence of 
rotation, suggesting a self-gravitating body. The appearance 
of the optical morphology is that of two gas rich disks col- 
liding almost perpendicularly, with the NIR image revealing 
the Heart nucleus to be likely associated with the horizon- 
tal disk component (the Wings) and the Body nucleus being 
more likely associated with the vertical structure comprising 
also the Tail - and potentially the Head. However, given the 
clearly discrepant velocity (and the significant stellar mass 
of the Head), it is difficult to imagine the Head being mate- 
rial merely ripped off a galaxy that originally had the Body 
component as its nucleus, and that component remaining 
such a conspicuous single nucleus. 

In fact, the optical images show a large obscuring dust 
region crossing perpendicularly the E-Wing, and also ex- 
tending some 6 kpc to the E and NE of the Head. This 
structure looks as if it could be separate from the two ma- 
jor horizontal and vertical components; and if this is true, it 
could be material from the same structure as the stellar com- 
ponent associated with the Head. The gas velocity (~ 15000 
km s _1 ) of this structure can be traced 4 kpc North out of 
the main Head area in NIR (PA=192; Fig [Till and together 
with the velocity at the eastern tip of the E-Wing (~ 14830 
km s _1 ), suggest a common kinematic origin with the Head 
(~ 14950 km s" 1 ). 

Though the final answer warrants detailed dynamical 
merger simulations out of the scope of this paper, we pro- 
pose that the data presented here strongly suggests a near- 
ULIRG merger of three pre-existing major components. The 
Bird would be, to our knowledge, the most-detailed studied 
case thus far. For a possible triple me rger case of a lower 
luminosty LIRG see lLfpari et al.l (|2000h ■ 
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4.3 AGN contribution and star formation rate 

We discuss next whether there is a possibility of an AGN 
contribution to the overall luminosity of the Bird, based on 
both optical and MIR considerations. First of all, there are 
no broad lines anywhere in the optical RSS spectra, nor do 
the line ratios suggest clear AGN activity. Some line ratios 
close to the HII vs. LINER dividing region come from ar- 
eas associated with outflows, and thus may be attributed to 
shock-heating. On the other hand, while most IRAC mid-IR 
colours of the Bird are those of normal spiral galaxies, along 
with strong PAH emission expected in starbursts contribut- 
ing to the 8 urn band, the Head is redder than the other 
components in all IRAC colours, suggesting a hotter dust 
contribution. 

To differentiate between AGN and star burst origin, 
we examine the IRS sp ec trum. Follow i ng Sturm et al.l 
|2002l ); IVerma et all (|2003h ; iFarrah et all (|2007h we calcu- 
late diagnostic line ratios: [NeIII]/[Ne II] = 0.11 ± 0.03, 
[SIV]/[SIII] = 0.6 ± 0.2 (the [SIV] line is totally within 
the wide silicate absorption feature), and [O IV] /[Nell] ~ 
0.03 ± 0.02 (the [OIV] line is a rather marginal detection). 
The Neon line ratio by itself, and also together with the Sul- 
phur line ratio, indicates a medium excitation starburst. The 
marginal detection of [O IV] raises the possibility of some 
AGN contribution, but the [OIV] /[Nell] ratio restricts this 
to definitely less than 5 per cent. In addition to the fine 
structure recombination line ratios, the strengths of the 6.2 
urn PAH feature and the 9.7 urn silicate absorption fea- 
ture (Ssil) place the Bird firmly in the s tarburst dominated 
ULI RG section of diagn ostic diagrams in ISooon et alJ (|2007l ) 
and lFarrah et al l l|2007l ). 

Thus, if there is any AGN contribution from the Head 
(and this would perhaps be the more unlikely location con- 
sidering the other two nuclei are more massive) it is negligi- 
ble compared to the total IR output of the Bird. 

We next estimate the star formation rate using 
the IR data available. First, we have fitted mod- 
els (lEfstathiou. Rowan-Robinson fc Siebenmorgenl |2000| ; 
lEfstathiou fc Rowan-RobinsonlbOOSt ) to the full IR SEP of 
the Bird, from NIR NACO and 2MASS data, to the FIR 
IRS, MIPS, and IRAS data. Fig. 1151 shows the data points 
and the SED model fit overplotted: the best-fitting model, a 
pure starburst, assumes an underlying exponentially decay- 
ing SFR with a burst in the last 7 x 10 7 years. The SFR at 
the peak of the burst was about 400 M©yr _1 or 188 MQyr -1 
when averaged over the last 7 x 10 7 years. 

Second, we utilize a newly calibrated relation between 
star formation rate and the PAH 6.2+11.2 urn luminosities 
Lp, based on correlations of thes e PAH features with total 
IR luminosity l|Farrah et al.ll2007f ): 



SFR [Moyr" 1 ] = 1.18 x 10" 



where Lp is the PAH 6.2+11.2 um luminosity in ergs 
s _1 . From the IRS spectrum (Fig. [5} we calculate Lp = 
1.6 x 10 43 ergs s _1 , i.e. we get a star formation rate of 
SFR = 192 M Q yr _1 , in very good agreement with the SED 
modelling. For refere nce, we calculate a rate of 120 M©yr _1 
from the traditional iKennicuttl (|l998h relation of SFR = 
1.72 x 10~ 10 Lir using the total IR luminosity of the Bird. 
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Figure 15. The 2MASS JHK, NACO K, Spitzer 3.6, 4.5, 5.6, 
8, 24, 70 um, and IRAS 12, 25, 60, 90 um photometric data 
overplotted with a model combining starburst and cirrus models. 
See text for details. 



4.4 Supernova rate 

Assuming core-collapse SN progenitors between 8 and 50 
Mq and a Salpeter initial mass f unction (IMF) with cut - 
offs at 0.1 and 125 Mq (see e.g. iMattila fc Meikld l200ll ). 
we can estimate the core-collapse SN rate of IRAS 19115- 
2124. Adopting the average SFR of 190 M yr _1 as ob- 
tained above, yields a core-collapse SN rate of about 1.3 
yr _1 . Comparison of our NACO K-b&nd images separated 
by 153 days did not reveal any new sources within the galaxy 
nuclear regions (see Section [3.1.3P down to limiting magni- 
tudes of K — 21.3 — 21.5. We compared these with template 
K-band light curves of ordinary (based on 11 type II and 
Ib/c SNe) and slowly-declining core-collapse SNe (based on 
the type HL SN 1 979C and the type Iln SN 1998S) from 
IMattila fc Meiklel (l200lj V_ Assuming the Galactic extinction 
law l|Rieke fc Lebofskilll985h we find that any slow decliner, 
suffering from an extinction of up to Ay ~ 30, should have 
been detected if exploding between the epochs of the two 
images or less than a year before the first epoch image. How- 
ever, in the case of ordinary core-collapse SNe and the same 
extinction, a detection would have been possible only near 
maximum light. If the extinction was considerably lower, Ay 
~ 5, we would have detected also ordinary events anytime 
between the epochs of the two images or up to 5 months 
before the first epoch image. 

Given the estimated core-collapse SN rate of 1.3 yr _1 
we estimate a Poissonian probability of 16 per cent for a 
SN non-detection assuming that all the core-collapse SNe 
within the nuclear regions of the Bird were slowly-declining 
and suffered from extinctions of less than Ay = 30. Recent 
observations of lu minous radio SNe wi t hin nuclear regions of 
nearby (U)LIRGs (fCoiina et alj|200ll: lLonsdale et al.ll2006l; 
l Alberdi et ail 120061 ; iParra et alJ 120071 ; iPerez-Torres et all 
[2007) show that these events have exploded within a dense 
circumstellar/interstellar environment. Therefore, a high 
rate of slowly-declining (in NIR) SNe might be expected 
in such environments. However, if most of the SNe within 
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LIRGs behave similarly to SNe observed outside nuclear 
starburst regions, the fact that no SNe were detected is not 
surprising, since we were sensitive to ordinary core-collapse 
SNe with high extinctions only near the epoch of maximum 
light. 

4.5 Evolution of the Bird as a merger 

We next discuss possible evolutionary sequences of the Bird 
system aided by numerical merger simulations. Detailed dy- 
namical modelling of a given observed galaxy merger is very 
time consuming as a large parameter-space of initial orien- 
tations and orbital configurations needs to be probed (e.g . 
Wahde fc Donnerll200ll;lNaab fc Burkerdl2003l ; lBarnes!l20o4 



Naab. Jesseit fc Burkerj 20061 ). Thus we here instead adopt 



an approach in which we model a typical 1:1 and a 3:1 galaxy 
merger in order to study the general predictions from merger 
simulations. 

The initial galaxy models are Milky Way-like galaxies 
with a disk mass of mo = 5.5 x 1O 1O M0 and a 20 per 
cent initial gas fraction, the bulge mass fraction is set to 
a third of the disk mass fraction and the galaxy is em- 
bedded in a dark matter halo with total virial mass of 
m v j r = 1.34 x 10 12 M(T) (see Johansson et al. 2007 in prep, 
ISpringel. Pi Matteo fc Hernquistl (|2005l ) for details) . For the 
3:1 merger the corresponding masses for the minor compo- 
nent are a third smaller than the values given above. The 
galaxies merge on a parabolic orbit with an initial separa- 
tion of r scp = 160kpc h" 1 and with a pericentric distance 
of r P eri = 5.0kpc h -1 ). The mergers are simulated with the 
paralell TreeSPH-code GADGET-2 <|SpringelU2005h and the 
simulations follow self-consistently the dark matter and gas 
dynamics, radiative gas cooling, star formation, as well black 
hole growth and the assoc iated feedback processes (Johans- 
son e t al. 2007 in prep., ISpringel. Di Matteo fc Hernquistl 
2005). Fig. [16] shows the resulting nuclear separations, rel- 
ative velocities and corresponding total relative star forma- 
tion rates of the two mergers as a function of time. 

Turning to our imaging data, the overall morphologi- 
cal appearance of the Bird, with multiple identifiable nuclei 
within a few kpc, and very obvious tidal tails, classifies it as 
a close 'pre-merger', i.e. a system after first approach and 
contact (when the tidal tails typically are generated), but 
before a full merger , when nuclei are fully coalesced (see e.g. 
IVeilleux et all 120021 ). As for the individual nuclei, the NIR 
light in the inner regions of the Heart show a surprisingly 
intact structure with bulge and disk components, and reg- 
ular appearing barred spiral arms extending from the inner 
300 pc out to some 2 kpc. The tidal Wings beyond this re- 
gion, however, break the regularity. Its rotation velocity vs. 
velocity dispersion ratio is in the range v ro t/<J ~ 2 — 3, i.e. a 
clear rotationally supported disk system. The nucleus of the 
Body, in contrast, shows a much more disturbed morphol- 
ogy, with a surface brightness profile intermediate to that 
of disk and spheroidal systems. There is clear evidence of 
disturbed spiral arm and/or bar structure within the cen- 
tral 1.3 kpc. It thus seems that the Body has thus far suf- 
fered more in the on-going merger, and looks to be already 
morphing towards a de Vaucouleurs stellar light distribu- 
tion. The v ro t/a ratio is ~ 0.4, which means that while it 
is a pressure supported system there is also significant ro- 
tational support involved. The ratio v ro t/& ~ 0.4 appears 



to be a typical value in (U )LIRG samples (e.g. IColina et al.l 
120051 : IDasvra et alj[2006ah . and the end result of the Bird 
merger will not be a slowly rotating massive elliptical, which 
are thou ght to be products of higher redshift galaxy forma- 
tion (e.g. lNaab et al.ll2007n . 

The relative (radial) velocity between the Heart and 
the Body was found to be very small, of the order of 50 
km s _1 , while the Head has a large positive velocity offset 
~400 km s" 1 from the other two. Since also the strongest 
star formation is from the Head (see next section), a pos- 
sible scenario is raised in which the Head is a foreground 
object falling into, or just passing, the Body/Heart system. 
Comparing to Fig. [16] we see that relative velocities in ex- 
cess of 400 km s _1 are possible for both 1:1 and 3:1 mergers; 
the Body/Heart:Head system would more closely resemble 
a 3:1 merger. The peak in the relative velocity is associated 
with the first passage of the galaxies, which also produces a 
clearly defined peak in the relative star formation rate at the 
same time of t ~ 1.3 Gyr (we did not attempt to model the 
level of SFR here, merely the pattern as a function of time). 
It is thus possible that the star formation of the whole sys- 
tem is dominated by the first triggering of star formation in 
the Head component, while the other nuclei are further into 
the merger and in between starburst phases. The locations 
of the bright NIR star clusters within the very extincted 
regions of the Head, at the edge closest to the other com- 
ponents, might be an indication of just such infall of the 
component through the tidal ISM. One of the weak points 
of this scenario is the strikingly normal appearance of the 
Heart nucleus. However, given that the Head is lower in mass 
by at least a factor of three than the Heart, it is plausible 
that the Head is experiencing stronger tidal effects compared 
to the Heart, leaving it more or less undisturbed. 

Using the masses for individual components from Sec- 
tion 13.41 the minimum mass for the end result, derived 
from dynamics, would correspond to a ~ 0.7 — 1.0 m* 
galaxy; total system mass from M/L relations is of the or- 
der ~ 2 . The end result masses are quite high, though 
not exceptional. Numerous studies (e.g. iGenzel et al.|l200ll 



IColina et aDl2005l : IDasvra et aT]|2006bl : iHinz fc Riek3l2(B 

indicate the end products of gas-rich spiral mergers, seen as 
LIRGs and ULIRGs, will be intermediate mass sub-m* el- 
lipticals. The IRAS 19115-2124 system is thus towards the 
high end of masses of local IR bright systems, but never- 
theless still well below the classification of a giant elliptical 
> 5 x lO n M . 

Finally, there are numerous bright concentrations 
within the Bird, and we have not studied them in de- 
tail in this work. We note, however, that many of them 
might be candidates for being, or becoming, tidal dwarf 
galaxies. As an example the blue concentration within the 
Tail, 3.5 arcsec South of the Body, is compact, distant, 
and massive (from NIR light) enough to satisfy the sur- 
vival criteria for tidal dwarfs (|Monreal-Ibero et al.l 120071 : 
IWetzstein. Naab. fc Burkertll2007h . 



4-5.1 Evolution of the Bird as a LIRG 

As far as the (U)LIRG phenomenon is concerned, the ac- 
cummulated wisdom indicates that the source of the huge 
IR luminosity comes from either AGN or nuclear star for- 
mation, or both, and in the case of the Bird we have shown 
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Figure 16. Simulations of a 1:1 merger (black curves) and a 3:1 merger (red curves). The left panel shows the nuclear separation, the 
middle panel the relative velocity of the galaxy nuclei and the right panel the relative star formation rate of the total system, as a 
function of time. The filled circle indicates the time of final coalescense of the two galaxy nuclei. 



that the star formation ovewhelmingly dominates. For both 
cases the idea is that the central engine is fed by infalling 
gas from the gas-rich progenitor d isks, and that this infall is 
greatly facilitated by b ars (see e.g. lMihos fc Hernquist 1996; 
|Pi Matteo et al.|[2007l . and references). 

Somewhat surprisingly, however, the majority of MIR 
output of the Bird system is not coming from either the 
major NIR bright galaxy nucleus (Body) with its embedded 
bar-like structure, or from the barred spiral nucleus (Heart). 
The IRAC and MIPS imaging clearly indicate that the Head 
is the source of the bulk of MIR flux, with secondary con- 
tribution from the Heart, while the contribution from the 
Body is negligible. 

The reason for non-conspicuous star formation in the 
major nuclei is not the dearth of fuel for star formation; the 
CO observations clearly show that there is an ample reser- 
voir of molecular hydrogen in all of the major components 
(Fig. I12|l . This would indicate that none of the nuclei have 
yet experienced the second, or later, star formation peaks 
expected from merging models (see Fig llGI) . 

There are recent reports of major off-nuclear and disk 
(rather than central ly concentrated) star formation in in- 
teracting galaxies, l|Lfpari et alj |2000| ; IWang et alj |2004| ; 
Ijarrett et al.l I2006I ). and even cases of non-enhanced or 
very mo derate star formation connected with interac- 
tions (e.g. lBergvall, Laurikainen fc Aaltoll2003l :l Cullen et al.l 
120071 ). These studies highlight the fact that the triggering 
mechanism of star formation in interactions is still far from 
understood. All these cases are, however, in lower luminos- 
ity systems. The Bird is thus quite unique in the sense that 
it is so close to being classified as an ULIRG, and still is not 
dominated by a central starburst in the major components. 
Or, if the Head component were to be considered as a major 
nuclear location, the starburst is not associated with those 
gas-rich locations with bar structures. 

That strong star formation is present in the Head need 
not, of course, be surprising in itself; as discussed above, 
we may be witnessing the Head interaction during its first 
passage by the two other nuclei that are further along their 
merging sequence. It is the timing which is more unique - 
ULIRGs are typically observed in later stages of interactions, 
and first passage starburst mom ents with associated large 
relative velocities are rarely seen l|Murphv et al.ll200ll ). 

In summary, the existence of large gas reservoirs and 
connected bar structures in merging disk galaxies do not 



automatically imply a central star-burst in that location. 
On the other hand, general merger modelling does predict 
epochs of lower SFR in between more intense episodes. In 
the case of the Bird it seems that the more massive nuclei 
are in such a phase, while the least massive component is 
caught virtually in the act of first high-speed fly-by. It would 
be very interesting to investigate further what are the fac- 
tors inhibiting the (U)LIRG phase in the Body nucleus, for 
example, or why exactly it has switched off. 



5 SUMMARY OF RESULTS AND 
CONCLUSIONS 

We have studied in detail the spectacular luminous infrared 
galaxy IRAS 19115-2124, which we have dubbed the Bird. 
Using near-diffraction limited adaptive optics imaging in the 
K-band (VLT/NACO), and matching space-based optical 
imaging (HST/ACS), combined with ground-based optical 
spectroscopy (SALT/RSS), we find the Bird to be a pre- 
merger of three components. 

Two of the components have unambiguously identifi- 
able galaxy nuclei. One has the quite regular composition of 
a barred spiral disk galaxy (Heart) while the other (Body) 
has a more disturbed nuclear region, with hints of both spi- 
ral/arm structure and a surface brightness profile interme- 
diate to that of disks and ellipticals. The third component 
(Head) appears to be an irregular galaxy - our kinemati- 
cal and dynamical data, however, strongly suggest that it is 
not merely a tidal dwarf borne out of the interaction of the 
other two, but rather a pre-existing less massive galaxy inti- 
mately involved in the interaction. The AO NIR imaging is 
crucial in classifying sources and nuclei due to the extremely 
clumpy nature of sources of this kind in the optical. In ad- 
dition, we see strong tidal tails (the Wings) and numerous 
concentrated knots, some of which may be candidate TDGs, 
and find a bright luminosity function for candidate SSCs. 

The combination of high-resolution imaging and spec- 
troscopy has resulted in the refinement of the systemic ve- 
locity of this LIRG, 14576±9 km s" 1 , and we have listed the 
kinematics of all the components. We have also investigated 
the line ratios, equivalent widths, and profiles of the impor- 
tant emission lines and find them to show characteristics of 
an H II type galaxy. 

We identify blueshifted Nal D absorption doublets 
mainly from the Heart, Head, and Body components. We 
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interpret these to be due to gas outflows from the Bird, with 
velocities typically in the range 100 - 300 km s _1 , but pos- 
sibly some flows ranging up to 600 km s -1 . Consistently, we 
also identify elevated [Nil] /Ha levels in blueshifted wings 
of emission lines from the nuclear components, suggesting 
shock-heating associated with outflows. 

By utilizing Spitzer IRAC and MIPS imaging, we are 
able to pinpoint where the MIR flux and the non-stellar 
emission is originating from. The Heart and especially the 
Head components are the regions producing most of the MIR 
flux. There is a hotter dust component affecting the 4 — 8[im 
bands at the location of the Head: Spitzer IRS spectroscopy 
is used to confirm the classification suggested by our op- 
tical spectroscopy, that this location, and the Bird system 
as a whole, is a star formation dominated system, with any 
AGN contribution below 5 per cent. The star formation is 
estimated to be ~ 190 Moyr -1 , corresponding to a core- 
collapse supernova rate of ~ 1.3 yr _1 . We detected no super- 
novae brighter than K ~ 21 in NIR images separated by 5 
months. 

We estimate averaged extinction over the system to be 
of the order Ay ~ 3 — 5, while the extinction derived from 
MIR diagnostics of the deeply embedded nuclear system re- 
sults in Ay « 30. Gas masses inferred are up to a level of 
M ~ 2 x 1O 9 M . 

We have determined the component masses of the Bird 
system using rotational as well as central velocity dispersion 
information. We note that without the high-resolution NIR 
imaging it would have been very difficult to differentiate be- 
tween velocity components related to nuclear bulges from 
those originating e.g. from HII regions in tidal regions or 
outflow components. Using conservative deblending of com- 
ponents we arrive at masses of 0.27 m,, 0.36 m* and 0.08 m* 
for the Body, Heart, and Head, respectively. Alternatively, 
using K-band light as a tracer, we find masses of 0.40 m*, 
0.25 m* and 0.16 m* for the inner 1 kpc regions of the same 
three regions, and 2.1 m* for the total system. 

In summary, most optical, NIR, MIR, and FIR emis- 
sion characteristics, as well as kinematical and dynamical 
attributes place IRAS 19115-2124 among typical ULIRG 
samples, even though it is not technically quite a ULIRG 
with an IR luminosity of Lm — 10 11 ' 9 L©. The main differ- 
ences are that it appears to be a three-component merger, 
and that its total dynamical mass is of the order of ~ lm*, 
or more, i.e. at the upper end of typical ULIRG merger rem- 
nants. Finally, we find the likely origin of the ULIRG phase 
engine to be the Head, the least massive of the three merger 
components, and the only one without signs of spiral arms 
and bars, which often are regarded as tell-tale signs of in- 
flowing gas to central starburst regions. A possible sequence 
of events is such that the Heart and Body met first; the Body 
especially is already in the process of morphing into an el- 
liptical like structure. The Head is observed during its first 
high-speed approach to the system, making it the strongest 
starburst at this time - the general characteristics of this 
scenario are consistent with numerical merger simulations. 
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